The Optimal Cosmic Epoch for Precision Cosmology 
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The statistical uncertainty in measuring the primordial density perturbations on a given comoving 
scale is dictated by the number of independent regions of that scale that are accessible to an observer. 
This number varies with cosmic time and diminishes per Hubble volume in the distant past or future 
of the standard cosmological model. We show that the best constraints on the initial power spectrum 
of linear density perturbations are accessible (e.g. through 21-cm intensity mapping) at redshifts 
2 ~ 10, and that the ability to constrain the cosmological initial conditions will deteriorate quickly 
in our cosmic future. 



Introduction. The early Universe was characterized 
by linear density perturbations with a fractional ampli- 
tude \8(t)\ <C 1, thought to have been seeded by quantum 
fluctuations during cosmic inflation across a vast range of 
scales spanning more than ~ 26 orders of magnitude |l| . 
As the Universe evolved, these perturbations were pro- 
cessed by its radiation and matter constituents. Recent 
surveys restricted to a small fraction of the total observ- 
able volume of the Universe allowed observers to read 
off cosmological parameters from the "Rosetta stone" of 
these density perturbations to an exquisite precision of a 
few percent [2|-|4j. 

In this paper we consider the fundamental limit to 
the precision of cosmological surveys as a function of 
cosmic time. Our analysis provides a global perspec- 
tive for optimizing future observations, and for assessing 
the ultimate detectability limits of weak features such as 
non-Gaussianity from inflation [5|. Existing cosmologi- 
cal data sets are far from optimal. For example, the pri- 
mary anisotropies of the cosmic microwave background 
(CMB) [J| sample only a two-dimensional (last scatter- 
ing) surface which represents a small fraction of the three- 
dimensional information content of the Universe at the 
redshift of hydrogen recombination, z ~ 10 3 . 

It is convenient to analyze the density perturbations 
in Fourier space with S^ = f d 3 r5(r) exp{ik • r}, where 
k = 2n/ A is the comoving wavenumber. The fractional 
uncertainty in the power spectrum of primordial density 

is given by [S, 01 , 

(i) 



perturbations P(k) = (|5k| 
AP(fc) 



P(k) ^NiF)' 

where the number of independent samples of Fourier 
modes with wavenumbers between k and k + dk in a 
spherical comoving survey volume V is, 

dN(k) = (2n)- 2 k 2 Vdk, (2) 

with N(k) being the integral of dN(k)/dk over the band 
of wavenumbers of interest around k. 

The maximum comoving wavelength A max that fits 
within the Hubble radius is set by the condition (see Fig- 
ure P, 

Amax(i) = 2i?H, (3) 



where i?H = c/(aH) is the comoving radius of the Hubble 
surface, a = (l + z) -1 is the scale factor corresponding to 
a redshift z (normalized to unity at the present time), and 
H(t) — (a /a) is the Hubble parameter at a cosmic time 
t. The corresponding minimum observable wavenumber 
fcmin(i) = 27r/A max naturally gives J mm dk[dN(k)/dk] s=s 
1. At z < 10 3 , H « H Q ^fl m a- 3 + fl A . Throughout the 
paper, we adopt a present-day Hubble parameter value 
of Hq = 70 km s _1 Mpc -1 (or equivalently h = 0.7) and 
density parameters f2 m = 0.3 in matter and fl\ — 0.7 in 
a cosmological constant [J]. 




FIG. 1: In the standard (post-inflation) cosmological model, 

a Fourier mode with a comoving wavelength A which enters 

the comoving scale of the Hubble radius _Rh = c(aH)^ 1 (in 

units of cHq 1 

to a redshift z 

radius at a later time (corresponding to a ox it). Hence, there is 

only a limited period of time when the mode can be observed. 



4.3 Gpc) at some early time (corresponding 
1), would eventually exit the Hubble 



Counting Modes. An observational survey is typically 
limited to a small fraction of the comoving Hubble vol- 
ume Vmax = -j-i?H, which evolves with cosmic time. Fig- 



ure[2]shows that V max starts small, then grows to a maxi- 
mum at the end of the matter dominated era, and finally 
diminishes due to the accelerated cosmic expansion. 
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FIG. 2: The Hubble volume Vmsxx., normalized by its present- 
day value V = ^-(c/Ho) 3 = 3.3 x 10 2 Gpc 3 , as a function of 
cosmic time t (top axis, in units of the Hubble time H^ 1 = 14 
Gyr) and scale factor a — (1 + z)" 1 (bottom axis). 



For the standard cosmological model in which the mat- 
ter density is dominated by cold dark matter (LCDM), 
nonlinear structure develops first on small spatial scales 
(large k) where it erases memory of the initial conditions. 
We associate the maximum wavenumber fc max for which 
the initial conditions are still in the linear regime with 
the minimum radius i?NL = ^mix °f a spherical top- 
hat window for which the root-mean-square amplitude 
of density perturbations is unity at the cosmic time of 
interest [8J|, 



<X 2 (i? NL ) = / 
JO 



4irk 2 dk 
(2tt) 3 



P(k) 



3ji(fci?NL) 



kR.NL 



= 1, (4) 



where ji(x) = (sinx — xcosa;)/a; 2 is the first spherical 
Bessel function. We adopt a present-day normalization of 
a(8h~ 1 ) = 0.8 for the standard LCDM power-spectrum 
with a spectral index n s = 1 [4|. 

Before nonlinear baryonic structure begins to form 
(z > 50) we associate the corresponding minimum wave- 
length, A m in = 27r/fc max , with the Jeans length for the 
baryons, Aj = 27r//cj, below which baryonic perturba- 
tions are smoothed out by gas pressure. We base this 
definition on the baryons since they are commonly used 
by observers to trace the underlying matter distribution. 
At z > 200 when the gas is thermally coupled to the 
CMB due to the residual fraction of free electrons after 



kj = 



2k B T lfi 

3/iTOp 



-1/2 



'fl m Ho = 0.35 kpc 



(5) 



where T 7j o = 2.73K is the present-day CMB tempera- 
ture, and /i = 1.22 is the mean atomic weight of neu- 
tral primordial gas in units of the proton mass m p . At 
z < 200 and before the first baryonic objects form, the 
gas temperature scales as (1 + z) 2 , and so kj scales up 
from the value in Eq. (J5J as [(l+z)/200] _1/2 . To keep our 
discussion simple, we ignore subtleties associated with 
the motion of the baryons relative to the dark matter |9| , 
the baryonic temperature fluctuations [lOj, and the X- 
ray and photo-ionization heating of intergalactic baryons 
by the first sources of light |ll| . These details (which are 
partially sensitive to the unknown properties of the first 
sources of light [8]) could reduce fc max at z > 10. Hence, 
our calculated & max should be regarded as the absolute 
upper limit, adequate for regions that are not influenced 
by these effects. 

Figure[3]shows the range of wavenumbers between fc m j n 
and fc max that are accessible within the Hubble radius as 
a function of cosmic time t and scale factor a(t). 
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FIG. 3: The range of comoving wavenumbers for which the 
linear power spectrum can be observed per Hubble volume as 
a function of cosmic time and scale factor. The minimum 
wavelength Amin = 27r/fc max is taken as the larger among 
the baryonic Jeans scale and the scale where nonlinear struc- 
ture forms at any given redshift. The maximum wavelength 
Amax = 27r/fc m i n is set by the Hubble diameter 2Ru- 

The maximum number of linear modes available to an 
observer at a cosmic time t is then given by, 



-Nmax(t) 



12tt 



(k 



3 

max 



k 6 



(6) 



The solid line of Figure [4] shows the resulting mini- 
mum fractional error in the power-spectrum amplitude 
A-P(fc max )/P(fc max ) = l/V^max, as a function of cos- 
mic time and scale factor. We have calculated Rnt, from 
Eq. flU using the LCDM power spectrum P(k) |J. The 
sharp rise in the minimum error at a > 0.1 occurs be- 
cause the nonlinear scale -Rnl increases rapidly above 
the Jeans length Aj around a redshift of z ~ 10. The rise 
is sharp since the root-mean-square amplitude of density 
fluctuations in LCDM <j(R) has a weak independence on 
small scales, implying that different small scales collapse 
at about the same cosmic time. Changing the definition 
of i?NL m Eq. @ to refer to 2-cr perturbations reaching 
an overdensity of unity [i.e., 2cr(i?NL) = 1] would simply 
shift this sharp rise to a value of a that is smaller by a 
factor of 2, since the linear growth factor of perturbations 
at z 3> 1 scales as a. 

A present-day observer located at a fixed vantage point 
can see multiple Hubble volumes of an earlier cosmic time 
t, which were causally disconnected from each other at 
that time. The total number of such regions available 
within a spherical shell of comoving width Ar = 2Rh (t) 
is .^regions = (47rr 2 Ar)/V max , where 



log(H t) 



r(t) 



cdt' 

wr 



(7) 



is the comoving distance to the shell center for obser- 
vations conducted at the present time to- The reduced 
statistical uncertainty of 1 / \/ (N rcg i ons x A max ) is shown 
by the dashed line in Figure |4] Under the most fa- 
vorable conditions, the probability distribution of den- 
sity perturbations can be compared to a Gaussian form, 
p(S)d5 = (27rCT 2 )~ 1/,2 exp{— S 2 /2a 2 }dS, to within a preci- 
sion of ~ 10 -9 , several orders of magnitude better than 
the level required for detecting non-Gaussianity from in- 
flation [5| • Measurements of the gravitational growth of 
P(k) for a particular k at multiple redshifts could test for 
modifications of general relativity or hidden constituents 
of the Universe to a similar level of precision. 

Conventional observational techniques, such as galaxy 
redshift surveys, Lyman-a forest spectra, or 21-cm in- 
tensity mapping [y, 1 12L 1 13j j , are currently limited by sys- 
tematic uncertainties (involving instrumental sensitivity 
and contaminating foregrounds) to levels that are well 
above the ultimate precision floor presented in Figure 2] 
But as advances in technology will break new ground for 



more precise measurements F4J , they might offer a qual- 
itatively new benefit - enabling observers to witness the 
evolution of cosmological quantities [such as P(k,z)] in 
real time over timescales of years [ljj |l_6| . 

Conclusions. Figure [4] implies that the most accurate 
statistical constraints on the primordial density pertur- 
bations are accessible at z ~ 10, when the age of the 
Universe was a few percent of its current value (i.e., hun- 
dreds of Myr after the Big Bang). The best tool for 




FIG. 4: The minimum fractional error attainable for the 
power-spectrum amplitude 1/\A/Vmax per Hubble volume, as 
a function of cosmic time and scale factor (solid line). The 
dashed line includes the reduction in the statistical uncer- 
tainty for a present-day observer who surveys a spherical shell 
of comoving width 2Rn(t) centered at the corresponding cos- 
mic times. 



tracing the matter distribution at this epoch involves in- 
tensity mapping of the 21-cm line of atomic hydrogen 
[g, [l2j, [l3| . Although the present time (a = 1) is still ad- 
equate for retrieving cosmological information with sub- 
percent precision, the prospects for precision cosmology 
will deteriorate considerably within a few Hubble times 
into the future. For simplicity, our quantitative results 
were derived for a Universe with a true cosmological con- 
stant. However, the ultimate loss of information holds for 
any type of accelerated expansion, even if the dark energy 
density is evolving in time. 

For pedagogical purposes, we considered the instanta- 
neous number of modes available on a space-like hyper- 
surfacc of a fixed cosmic time t, under the assumption 
that a typical cosmological survey would focus on a small 
fraction of V max (t) in which the constant time approxi- 
mation is valid. Otherwise, the evolution of the density 
field and its tracers needs to be taken into account. 

Since the past lightcone of any observer covers vol- 
umes of the Universe at earlier cosmic times, one might 
naively assume that the accessible information only in- 
creases for late-time observers, as past information will 
stay recorded near the horizon even in the distant future 
13, |23[ [24 1 . However, in practice the exponential expan- 
sion will erase all this information in the future [17H22|. 
Beyond a hundred Hubble times (a trillion years from 
now) , the wavelength of the CMB and other extragalac- 
tic photons will be stretched by a factor of > 10 29 and 



exceed the scale of the horizon [25| (with each photon 
asymptoting towards uniform electric and magnetic fields 
across the Hubble radius), making current cosmological 
sources ultimately unobservable. While the amount of 
information available now from observations of our cos- 
mological past at z ~ 10 is limited by systematic uncer- 
tainties that could potentially be circumvented through 
technological advances, the loss of information in our fu- 
ture is unavoidable as long as cosmic acceleration will 
persist. 
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